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Introduction 

Two  problems  are  encountered  in  designing  nozzles  for  super-sonic  l2  injection  in  the 
chemical  oxygen  iodine  laser  (COIL).  First,  prediction  of  the  mixing  dynamics  using 
computational  fluid  dynamics  (CFD)  methods  is  a  challenging  task.  Madden  has  shown  that 
high-level  CFD  models  are  needed  to  capture  the  complex  turbulent  flow  characteristics. 
Secondly,  uncertainties  in  the  I2  dissociation  mechanism  make  prediction  of  the  dissociation  rate 
and  the  concomitant  energy  cost  of  dissociation  difficult  to  assess.  Validation  of  the  CFD  models 
for  super-sonic  injection  is  a  critical  task  in  the  sequence  that  leads  to  the  development  of 
efficient  super-sonic  mixing  nozzles.  The  group  at  AFRL  (Kirtland  AFB)  is  using  planar  laser 
induced  fluorescence  (PL  IF)  of  the  I2  B-X  transition  to  image  the  flow  in  a  Mach  2.5  nozzle  that 
is  equipped  with  super-sonic  injectors.  The  images  observed  in  their  current  experiments  show 
turbulent  structures  and  poor  penetration  of  the  secondary  flow  (attributed  to  the  shallow  angle  of 
injection).  The  raw  PLIF  images  provide  a  qualitative  picture  of  the  1 2  spatial  distribution. 
However,  the  l2  excitation  and  relaxation  kinetics  must  be  taken  into  account  in  order  to  convert 
the  images  into  quantitative  concentration  maps.  The  objective  of  the  program  described  here 
was  to  investigate  the  quenching  and  energy  transfer  kinetics  of  I2(B)  in  a  supersonic  flow,  and 
to  develop  kinetic  models  that  can  be  used  in  the  analysis  of  PLIF  data. 
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Quenching  cross-section  measurements 

a.  Characterization  of  the  nozzle 

Quenching  of  hfB)  by  N2,  O2  and  He  over  the  temperature  range  from  295  to  1 50  K  was 
investigated  in  a  Mach  2.6  nozzle.  The  first  task  of  this  project  was  to  configure  the  nozzle  and 
verify  that  supersonic  flow  conditions  had  been  achieved.  Fig.  1  shows  the  nozzle  assembly 


Fig  1 .  Photograph  of  nozzle  assembly. 


As  can  be  seen  in  this  picture,  the  nozzle  was  equipped  with  several  ports  to  permit  local 
pressure  measurements.  The  port  on  the  leftmost  side  of  the  picture,  tapped  into  the  bottom  of 
the  assembly,  was  used  to  measure  the  pressure  in  the  plenum.  The  first  port  (leftmost)  on  the 
top  of  the  nozzle  accessed  the  transonic  region  of  the  flow.  Eight  other  ports  sampled  the 
supersonic  region  of  the  flow.  Lastly,  the  pressure  tap  mounted  in  the  transparent  side-wall  of 
the  nozzle  assembly  provided  the  Pitot  pressure  at  the  end  of  the  expansion.  The  critical 
dimensions  of  the  nozzle  assembly  are  shown  in  Fig  2. 
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Fig.  2.  Cross-section  of  the  nozzle  assembly.  The  pressure  could  be  measured  in  the  subsonic 
region  (Pi),  and  in  the  transonic  region  (P2)  and  at  intervals  of  I  cm  along  the  flow  (P3-P10).  A 
Pitot  tube  was  used  to  measure  total  pressure  at  the  nozzle  exit. 


Supersonic  operation  of  the  nozzle  was  verified  for  flows  of  air,  and  the  pure  gases  N2,  O2  and 
He.  A  typical  pressure  profile  for  the  nozzle  operating  with  a  mixture  of  h  (0.3  %)  in  N2  is 
shown  in  Table  I  below 


Table  1 .  Pressure  distribution  along  the  flow  for  N2/I2  mixture 
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This  profile  was  consistent  with  supersonic  gas  flow  with  a  terminal  Mach  number  of  2.66. 

Local  temperatures  within  the  jet  were  determined  using  both  the  thermodynamic 
relationship 
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(where  y  is  the  heat  capacity  ratio)  and  by  fitting  to  the  rotational  line  intensity  distribution  of  the 
I2  B-X  spectrum  recorded  using  laser  excitation.  Fig.  3  shows  the  experiment  configured  for 
laser  induced  fluorescence  measurements. 


Fig.  3.  Schematic  view  of  the  LIF  measurement  system. 


Note  that  the  small  bypass  cell  shown  on  the  right  hand  side  of  this  figure  was  used  to  measure 
the  concentration  of  I2  in  the  flow  using  the  B-X  absorption  band  near  500  nm.  The  dye  laser 
beam  was  directed  along  an  axis  perpendicular  to  the  gas  flow.  Laser-excited  fluorescence  was 
detected  by  looking  along  the  flow  axis.  The  laser  system  used  for  these  experiments  consisted 
of  a  Nd/YAG  pumped  dye  laser,  operated  at  wavelengths  near  577  nm.  The  linewidth  of  this 
laser.  0.06  cm'1,  was  sufficiently  narrow  to  permit  rotational  resolution  of  the  B-X  bands. 
Examples  of  rotationally  resolved  spectra  recorded  in  the  nozzle  are  shown  in  Fig.  4.  I  lere  the 
red  trace  was  recorded  with  the  laser  positioned  in  the  plenum,  were  the  local  temperature  was 
that  of  the  laboratory.  The  green  trace  was  recorded  with  the  laser  positioned  6  cm  downstream 
from  the  throat  of  the  nozzle.  Rotational  temperatures  were  derived  from  these  data  using  a 
spectral  simulation  program  (c.f.  Van  Marter  et  al.  Chem.  Phys.  Lett.  260,  201  (1996)).  Overall, 
the  agreement  between  the  temperatures  obtained  from  Eq.  (1)  and  those  obtained  from 
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Fig.  4.  LIF  spectra  obtained  using  N2/I2  with  P|=  50TorT,  Pe=  3.7  Torr,  Ppj,ot=32  Torr  and 

M=2.52. 

the  spectra  was  satisfactory  (errors  of  no  more  than  10%).  In  several  experiments  we  varied  the 
conditions  in  the  nozzle  by  adding  a  flow  of  air  to  the  region  downstream  of  the  nozzle  exit  plane 
(which  has  the  effect  of  decreasing  the  pumping  speed).  With  this  approach  we  could  hold  the 
pressure  constant  at  a  particular  point  in  the  flow  and  vary  the  local  temperature  by  varying  the 
plenum  pressure  and  the  pumping  speed.  To  verify  the  proper  operation  of  the  nozzle  under 
these  conditions  we  compared  the  Mach  number  derived  from  Eq.  I  (denoted  by  Madiab)  with  that 
obtained  using  the  pressure  measured  by  the  Pitot  tube  (Mpllol).  The  latter  was  obtained  from  the 
expression 


which,  for  y=l  .4,  reduces  to 
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Table  2  shows  a  comparison  of  the  Mach  number  at  P?  determined  from  Eq.’s  1  and  3.  Here  it 
can  be  seen  that  the  agreement  is  good  for  high  values  of  the  Mach  number  (>2)  but  there  are 
discrepancies  for  M  in  the  range  of  1.25-1.8.  For  this  range  we  found  that  the  Mach  numbers 
derived  from  the  Pitot  tube  pressure  were  more  reliable. 


Table  2.  Mach  numbers  determined  from  Eq.’s  1  and  3  for  P9=2  Torr 


Pi,  Torr 

Ppitot,  Torr 

Mplt0( 

Madiab 

43 

18.8 

2.63 

2.65 

37.6 

18.5 

2.62 

2.56 

32.8 

17.1 

2.51 

2.48 

27.4 

15.6 

2.4 

2.36 

22.5 

12.95 

2.16 

2.24 

20.2 

12.7 

2.15 

2.17 

16 

9.4 

1.82 

2.02 

12.6 

6.75 

1.5 

1.86 

10.6 

5.2 

1.25 

1.75 

(b)  Measurement  of  quenching  rate  constants 

Time-resolved  fluorescence  decay  measurements  were  used  to  characterize  the  12(B) 
quenching  rate  constants  as  a  function  of  temperature.  Fig.  5  shows  typical  time-resolved  signals 
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Fig.  5.  l2(B,v=14)  fluorescence  signals  for  I2/N2  Pt=2.5  Torr 


recorded  for  h  in  N2  with  the  probe  laser  positioned  on  the  flow  centerline,  directly  below  the  P7 
pressure  tap.  The  quenching  rate  constant  at  each  temperature  was  determined  from  the 
expression 

r,^  =  r,P+kM  (4> 

where  fllkt  is  the  measured  decay  rate,  r  is  the  spontaneous  decay  rate,  kq  is  the  quenching 

rate  constant,  and  [q\  is  the  number  density  for  the  quenching  gas  q.  Fig.  6  shows  a  plot  of  the 
rate  constant  for  quenching  of  l2(B,  v=14)  by  N2  over  the  temperature  range  from  1 10-295  K. 


Fig.  6.  Rate  constant  k(T)  for  quenching  of  hfB,  v=14)  by  N2 

Over  the  range  from  1 50  -295  K  the  rate  constant  is  approximately  linearly  proportional  to 
temperature,  which  indicates  that  the  cross  section  for  the  quenching  process  increases  slowly 
with  decreasing  temperature.  In  the  region  below  150  K  the  rate  constant  levels  off  at  an 
approximately  constant  value.  If  this  effect  is  real,  it  would  indicate  that  the  cross  section  was 
increasing  more  rapidly  with  decreasing  temperature.  This  could  be  a  consequence  of  the 
increasing  importance  long-range  attractive  forces  at  low  temperatures.  However,  there  was  also 
some  concern  that  the  lowest  temperature  measurements  were  subject  to  systematic  errors.  For 
temperatures  above  1 50  K  the  rate  constant  was  we 1 1 -represented  by  the  expression  (cm1  s'1) 
kNf(T j  =  fl.97±O.I9jxlO"n  +fl.59±0.09jxl0‘13  x  T 
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The  temperature  dependence  of  the  rate  constant  for  quenching  of  I2(B,v=I4)  by  O2  is 
shown  in  Fig.  7.  Here  again  it  can  be  seen  that  the  rate  constant  has  a  near  linear  dependence  on 
temperature  for  T>  1 50  K  and,  in  this  instance,  the  rate  constant  appears  to  increase  with 


Fig.  7.  Rate  constant  k(T)  for  quenching  of  l2{B,  v=!4)  by  O2 

decreasing  temperature  below  140  K.  For  T>1 50  K  the  rate  constant  is  given  by 
Ao/r;  =  f2.18  +  036;xl0-u+fl.89±0.19;xl0-,3x7’ 

Quenching  of  12(B)  by  !  le  was  found  to  be  very  slow.  Fig.  8  shows  a  decay  curve  recorded  at  a 
total  He  pressure  of  3  Tort  at  T~I60  K.  This  fluorescence  decay  rate  was  close  to  that  expected 
if  the  deactivation  was  determined  by  self  quenching  alone.  Few  measurement  were  done  with 
pure  He  as  the  pumping  speed  was  not  sufficient  to  achieve  stable  supersonic  operation  of  the 
nozzle.  Since  these  measurements  were  carried  out,  new  nozzle  blades  have  been  fabricated  that 
define  a  contour  that  has  a  small  throat  size,  and  this  will  reduce  the  gas  flow  with  He  to  a  level 
that  the  current  pumping  system  is  able  to  cope  with.  Further  measurements  of  the  quenching 
rate  constants  using  the  small  nozzle  will  be  carried  out  in  the  next  round  of  measurements. 
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Fig.  8.  Fluorescence  decay  12(B)  in  a  He  supersonic  flow  at  P=3  Torn 


Appendix.  A1AA  conference  paper  describing  the  initial  kinetic  model  developed  for  analysis  of 
[2  PLIF  images. 
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Introduction. 

Imaging  of  fluid  dynamic  flowfields  using  non-intrusive  diagnostics  js  a  matter  of  practical  interest  to  both 
experimentalists  and  theoreticians.  All  flow  diagnostics  such  as  Pitot  tubes,  hot-wire  probes,  pressure  taps,  panicle 
image  velocimeiry,  etc,,  offer  the  ability  to  obtain  valuable  information  about  the  flow  characteristics  but  perturb  the 
local  flow  stale  to  some  degree.  Evaluations  must  be  made  comparing  the  diagnostic  need,  the  type  of  information 
provided,  and  the  cost  of  the  diagnostic  among  other  considerations  to  determine  the  best  choice.  Non-intrusive 
diagnostics,  particularly  those  using  coherent  radiation  to  generate  a  secondary  radiative  response  from  the  particles  in 
the  gas,  are  attractive  because  of  the  inherent  ability  provide  local  measurements  with  minimal  perturbation  of  the  flow. 
In  the  case  of  chemical  oxygen-iodine  laser  (COIL)  flowfields,  laser  induced  fluorescence  (LIF)  techniques  are  attractive 
for  a  variety  of  reasons,  the  non-intrusive  nature  of  LIF,  the  ability  to  leverage  existing  COIL  development 
infrastructure,  the  ability  to  get  detailed  localized  information  about  flow  structure,  and  the  natural  combination  of  LIF 
with  the  presence  of  l2  in  COIL  flowfields,  termed  laser  induced  iodine  fluorescence  (LIIF),  A  drawback  to  the  use  of 
LIIF  in  COIL  flowfields  is  the  compressible  nature  of  the  flows,  with  large  changes  in  pressure,  temperature,  and  local 
species  densities  occurring  within  the  field  that  must  be  accounted  for  in  making  quantitative  interpretations  of  LIF 
generated  imaging.  The  development  of  a  model  for  the  12  fluorescence  that  accounts  for  the  thermodynamic  variations 
while  describing  the  spatial  variation  of  photon  production  is  useful  from  this  standpoint  as  it  allows  for  quantification  of 
the  compressibility  effects  from  an  interpretive  standpoint  while  at  the  same  time  providing  a  tool  for  validation  and 
verification  of  computational  fluid  dynamic  (CFD)  models. 

LIF  as  a  flow  imaging  technique  originates  in  the  work  of  Daily1  Wrobel,1  and  Levy3  who  applied  the  non-intrusive 
nature  of  the  technique  as  a  convenient  mechanism  for  imaging  the  development  of  flow  structure  for  a  variety  of 
disparate  Hows  including  mixing  layers,  reacting  flows,  and  expanding  jets.  Original  work  by  Rapagnani  and  Davis'* 
applied  the  technique  as  a  flow  diagnostic  for  the  interrogation  of  the  mixing  systems  in  chemical  lasers  as  a  steady  state 
flow  diagnostic  utilizing  a  cw  laser  source  exciting  the  X-^B  transitions  of  I2,  although  the  issues  associated  with 
thermodynamic  and  quenching  rate  constant  variation  within  the  flow  due  to  compressibility  were  not  addressed, 
McDaniel5’6  addressed  the  issue  of  use  of  LIIF  as  a  quantitative,  steady  state  diagnostic  for  compressible  flows.  His 
work  found  particular  application  in  the  mixing  systems  of  high  speed  engines  where  the  ability  measure  local  quantities 
over  small  volumes  was  particularly  useful  for  diagnosing  mixing  rates  in  these  systems  and  creating  multidimensional 
datasets  for  validation  and  verification  of  CFD  models.  The  issue  of  LIF  intensity  variation  resulting  from 
compressibility  effects  was  addressed  by  Fletcher7  and  Hartfield,*  again  for  cw  laser  sources  and  steady  state  diagnosis. 
Given  the  analysis  performed  by  McDaniel  in  developing  the  LIF  for  compressible  flow  application,  the  use  of  this 
model  on  CFD  datasets  was  an  attractive  extension  of  the  theoretical  development.  Hartfield  performed  such  work, 
applying  the  model  as  a  computational  fluid  imaging  technique  to  steady  state,  3-D  CFD  results  and  comparing  directly 
to  LIF  imaging  of  high  speed  mixing  system  experiments.  As  mentioned  earlier,  LIF  is  an  obvious  candidate  for 
application  to  COIL  flowfields  due  to  the  use  ofl2  as  a  reactant  in  these  devices.  Rapagnani  and  Davis9  re-applied  the 
technique  for  COIL  flowfield  imaging  as  a  validation  and  verification  tool  complementing  the  3-D  CFD  simulation 
efforts  by  Rapagnani.10  Recent  work  by  Noren 11,12  utilizing  LIF  of  I2  for  the  investigation  of  COIL  mixing  nozzles  is 
coupled  directly  to  this  effort,  and  within  this  context  the  current  effort  represents  a  validation  and  verification  study  for 
the  computational  component. 

The  work  performed  in  this  investigation,  while  focusing  on  LIFt  differs  from  the  work  of  McDaniel,  Fletcher,  and 
I  lartield  in  the  use  of  short  pulse  laser  pumping  to  excite  the  \2  as  opposed  to  the  cw  pumping  sources  utilized  in  that 
work.  Thus  the  limiting  assumptions  and  simplifications  that  result  from  the  use  of  a  cw  laser  source  no  longer  apply, 
and  the  ordinary  differential  equations  that  describe  the  production  and  loss  of  the  X  and  B  states  must  be  integrated 
directly  over  the  period  of  the  laser  pulse  to  track  photon  production  during  the  B^X  excitation  and  reemission  process. 
Additional  differences  regarding  the  temperature  dependence  of  the  quenching  processes  as  well  as  the  mechanisms  that 
determine  the  populations  of  the  vibrational -rotational  states  pumped  by  the  laser  pulse  within  the  ground  state  manifold 
will  be  elucidated  below. 

Problem  and  Methodology. 

The  experiment  that  is  modeled  here  is  the  planar  LIIF,  most  commonly  referred  to  as  PLIF,  interrogation  of  a  Mach 
3  A  He/12  jet  injected  into  a  He/02  flow  that  is  expanding  to  Mach  2,2.  The  point  of  injection  is  just  downstream  of  the 
primary  flow  nozzle  throat  on  the  upper  surface.  Figure  L  illustrates  the  experiment  configuration  consisting  of  the 
nozzle  flow  hardware,  the  laser  sheet  passing  through  the  gas  flow,  and  the  CCD  camera  and  Fig.  2  shows  the  secondary 
flow  injector  orifice  as  it  was  machined  into  the  primary  flow  nozzle  wall.  The  primary  flow  consists  of  505  mmol/s  He 
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and  125  mmol/s  02  with  a  total  pressure  of  65  torr  and  a  total  temperature  of  298.15  K.  The  secondary  flow  consists  of 
125  mmol/s  He  and  1  mmol/s  l2  with  a  total  temperature  of  425  K  and  a  total  pressure  of  453  torr. 

The  PL1F  pumping  laser  consists  of  a  tunable  dye  laser  using  Rhodamine  6G  dye  pumped  by  a  Nd:YAG.  The 
measured  wavelength  of  the  tuned  dye  laser  is  565.0947  nm  and  the  measured  bandwidth  is  22.53  GHz.  The  pulse  width 
for  the  dye  laser  output  is  30  ns.  The  output  beam  of  the  dye  laser  is  optically  shaped  to  form  a  600  fm  thick  sheet  prior 
to  transmission  through  the  experiment  for  fluorescing  the  I2.  The  light  collection  system  uses  a  photo- intensified 
5 1 2x5 1 2  pixel  array  CCD  camera  gated  to  20  ns. 

Relative  to  Fig.  I  the  laser  light  can  interrogate  the  gas  at  either  0°  or  45°  angles  when  the  angle  is  measure  with 
respect  to  an  axis  passing  through  the  experiment  parallel  to  the  primary  flow  direction.  The  0°  interrogation  is 
accomplished  by  passing  the  beam  upstream  through  the  gas  parallel  to  the  primary  flow  direction  through  a  Lexan 
window  downstream  of  the  test  section.  The  45°  interrogation  is  illustrated  in  Fig.  I ,  with  the  light  passing  upstream  and 
across  the  flowfield  with  the  45°  angle  relative  to  the  primary  flow  axis.  The  camera  is  positioned  orthogonal  to  the 
laser  light  sheet  for  photon  collection. 

The  approach  taken  for  the  development  of  a  model  for  the  PUF  imaging  of  the  flowfletd  within  the  experiment  is 
two-fold.  The  first  part  deals  with  the  development  and  execution  of  a  time  dependent,  3-D,  Navier-Stokes  simulation 
for  a  LI  IF  experiment.  The  second  part  deals  with  the  development  and  execution  of  a  model  for  the  physical  processes 
underlying  LI  IF  in  this  experiment.  What  is  implicit  here  is  the  de-coupling  of  the  physical  model  for  the  fluid  dynamics 
from  the  physical  mode)  for  LIIF.  This  approach  is  justified  based  upon  the  difference  in  time-scales  for  the  pulsed  UIF 
and  the  fluid  dynamics.  The  peak  velocity  in  the  primary  flow  direction  within  the  experiment  is  order  1600  m/s.  Over 
the  course  of  the  period  during  which  photons  are  collected  by  the  CCD,  the  camera  gate  of  20  ns,  the  flow  in  this  peak 
velocity  region  will  translate  3.2-  10  s  m.  Given  the  length  scale  of  the  exit  diameter  of  the  jet  of  order  0.001  m,  a 
representative  length  scale  for  the  injected  jet  fluid,  the  flow  translation  distance  of  3.2  I0'S  m  represents  3.2%  of  this 
value.  Practically  speaking,  the  flow  is  frozen  during  the  period  of  the  camera  gate  and  for  this  reason  the  model  for  the 
PL  IF  processes  can  be  de-coupled  from  the  CFD  results. 

The  CFD  model  for  the  flow  experiment  utilizes  a  numerical  solution  of  the  Navier-Stokes  continuity  equations  for 
mass,  momentum,  and  energy,  with  individual  mass  continuity  equations  for  each  of  the  species.  It  is  assumed  that  no 
chemical  reactions  occur  within  the  gas  during  the  experiment  and  formation  of  I2-He,  I2-02l  or  Ij-lj  van  der  Waals 
complexes  is  assumed  negligible  as  well.  Accurate  description  of  the  molecular  diffusion  of  the  species  within  the  gas  is 
important  for  the  compressible,  low  Reynolds  number  conditions  of  these  flows  especially  when  considering  the  desire 
to  compare  the  mixing  of  the  secondary  fluid  as  indicated  by  the  PLIF  model  imaging  of  the  CFD  simulation  to  the  PL1F 
images  from  the  experiment.  The  Ramshaw-Dukowicz13  approximation  for  multi -component  diffusive  transport  is  used 
including  both  gradient  and  pressure  diffusion  terms.  The  physical  domain  of  the  primary  and  secondary  flows 
represented  in  Figs.  1  and  2  is  represented  in  the  computational  domain  of  the  simulation  using  a  13  block,  51.2  million 
cell  grid  shown  in  Fig.  3.  The  grid  compression  at  the  no-slip  surfaces  uses  a  minimum  spacing  at  the  wait  set  to  y  -  1. 
The  primary  and  secondary  subsonic  inflow  boundaries  for  the  simulation  fix  the  total  pressure  and  temperature  white 
setting  the  second  derivative  of  the  static  pressure  to  0,  allowing  pressure  waves  to  pass  through  while  fixing  the  mass 
flux.  The  outflow  boundary  condition  is  a  supersonic  outflow  and  sets  the  second  derivative  of  the  dependent  variables 
to  zero.  The  solid  surfaces  are  treated  as  no-s!ip  boundaries  with  the  temperatures  of  the  lower  nozzle  surface  and 
sidewalls  set  to  a  temperature  of  298. 1 5  K  and  the  heated  upper  surface  in  which  the  He/I2  orifice  is  located  set  to  400  K. 

A  consequence  of  the  frozen  flow  approximation  is  that  the  inviscid  and  molecular  transport  aspects  of  the 
continuity  relationships  described  by  the  Navier-Stokes  equations  become  negligible  and  the  only  terms  that  remain 
significant  are  the  time  derivatives  of  the  various  upper  and  lower  states  in  the  I2  X  and  B  vibrational-rotational 
manifolds  and  the  source  terms  that  describe  their  production  and  loss.  The  physical  processes  that  occur  within  the 
fluorescing  process  of  I2  can  be  organized  as  follows:  the  interaction  of  the  laser  light  with  the  ground  electronic  state  of 
l2  such  that  a  specific  vibrational-rotational  state  within  the  total  manifold  of  X  states  is  depopulated  and  a 

specific  vibrational -rotational  state  ( v',J' )  within  the  total  manifold  of  the  l2  B  electronic  state  is  populated;  through 
direct  collisions  and  long  range  interactions  with  other  molecules  the  ground  v",J"  state  is  re-populated  and  the  v',J' 
excited  state  is  de-populated;  through  spontaneous  and  simulated  emission  the  entire  manifold  of  populated  B  states  is 
de-populated.  Mathematically  this  is  described  as  the  set  of  coupled  ordinary  differential  equations: 
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In  this  $ei  of  equations,  a**,  is  the  absorption  cross-section  for  the  transition  being  pumped  and  p(v)  is  the  photon  flux 
per  unit  area.  gx-  and  gB ■  are  the  degeneracies  for  the  ground  and  excited  states  respectively.  X'  is  the  population  of  the 
ground  vibrational -rotational  state  that  is  pumped  to  the  excited  vibrational -rotational  state  B*  within  the  B  electronic 
state  manifold.  X  is  the  sum  total  population  of  all  the  remaining  vibrational -rotational  stales  within  the  ground 
electronic  state  manifold  and  B  is  the  sum  total  population  of  all  the  remaining  vibrational -rotational  states  within  the  B 
electronic  state  manifold,  kmX  is  the  total  rate  at  which  the  X  states  are  depopulated  to  fill  X*  during  the  pumping 
process  through  interactions  with  the  He  and  02  molecules  in  the  gas  and  as  such  is  a  rotational  energy  transfer  rate 
constant.  Similarly,  k,nB  is  the  total  rate  constant  for  the  B  states  are  colJisionally  depopulated  to  fill  B\  Along  these 
same  lines,  kouf  X  is  the  rate  constant  for  population  transfer  from  X'  to  the  X  states  and  k^g  is  the  rate  constant  for 
transfer  from  the  Br  state  to  the  B  states.  F^  is  the  radiative  decay  rate  of  the  B  (and  B ')  state,  Fdiss  is  the  rate  at  which 
the  B  state  levels  are  lost  through  spontaneous  (electronic  p  re  dissociation)  and  collision  induced  dissociation  processes. 
For  the  conditions  of  the  experiments  described  here,  F^  is  given  by 


Cfa*  =  rtP  +  kqH*He  +  k 02  +  k ^  /  2 


m 


where  the  kqM  terms  are  the  rate  constants  for  quenching  by  collision  induced  dissociation.  rsp  is  the  decay  rate  for 
spontaneous  predissociation.  This  results  from  the  coupling  of  B  state  levels  with  a  dissociative  state  that  correlates  with 
ground  state  atomic  products.  The  states  are  coupled  by  rotational  motion  of  the  molecule,  such  that  the  predissociaiion 
rate  is  proportional  to  the  rotational  energy.  The  total  decay  rate  for  levels  of  the  B  state  is  then  given  by 
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Values  for  the  radiative  and  predissociative  decay  lifetimes  (reciprocals  of  the  decay  rates)  used  here  are  shown  in  Table 
2.  In  addition,  Table  2  also  lists  the  collision-free  lifetimes  (to),  which  reflect  the  combination  of  radiative  and 
predissociative  decay* 

Implicit  within  this  model  is  the  functional  dependence  on  pressure  and  temperature,  as  many  of  the  terms  in  these 
equations  have  some  dependence  on  the  two  variables.  The  rate  constants  for  rotational  energy  transfer,  Jt„  include 
temperature  dependence  and  the  product  of  ktM  includes  both  temperature  and  pressure  through  the  number  density.  The 
12  dissociation  rate  is  scaled  by  the  temperature  through  the  collision  frequency  of  l2  with  itself,  He,  and  02.  The 
absorption  cross  section  <rahs  is  a  weak  function  of  pressure  and  temperature  and  the  initial  population  of  the  ground 
vibrational-rotational  state  is  a  function  of  temperature  through  the  Boltzmann  population  distribution.  The  temperature 
variation  of  the  distributions  for  the  odd-./  levels  within  v=G  of  the  ground  electronic  state  of  l2  is  shown  in  Fig,  4, 
illustrating  the  wide  variation  of  the  equilibrium  population  distribution  within  the  range  of  temperatures  expected  m  the 
flowfield  simulated  here. 

The  temperature  dependencies  of  the  l2(£)  quenching  rate  constants  were  examined  by  Azyazov  and  Heaven14  in 
nozzle  flow  experiments  using  He  and  02  as  the  buffer  gas.  The  use  of  the  nozzle  flow  permitted  temperature  variation 
studies,  and  measurements  were  made  between  150  K  and  295  K.  Over  this  range  the  rate  constants  exhibited  a  7^  ; 
dependence.  The  results  were  consistent  with  Cappelle  and  Broida’s18  expression  for  the  collisional  part  of  F^  in  terms 
of  temperature  independent  cross  sections: 
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with  //  being  the  reduced  mass  for  the  collisions  with  lz  and  A  being  the  Boltzmann  constant.  The  collisions]  cross 
section  for  self-collisions*  <7^,  is  70'  I O’20  m\  for  collisions  with  He  o^=0,30- 1020  m2,  and  for  02  cr^  =5  10  20  m3. 

Knowledge  of  the  degeneracy  of  the  individual  states  within  the  iodine  molecule  X  and  B  manifolds  was  required  for 
calculations  of  excitation  rates  and  population  distributions.  For  any  rotational  level  of  a  diatomic  rotor  the  degeneracy  is 
(2 J+l)  for  the  rotational  angular  momentum.  Further  degeneracy  for  I2  is  derived  from  the  nuclear  spin  states  of  the 
molecule.  The  nuclear  spin  of  127I  is  5/2*  so  that  the  total  spin  quantum  number  for  the  molecule  takes  the  values  0f  !t  2, 
3*  4*  and  5,  Due  to  the  restrictions  imposed  by  nuclear  exchange  symmetry,  the  even  nuclear  spin  states  combine  with 
even -7  states  while  odd-spin  states  combine  with  odd-7  for  the  ground  state*  which  has  gerade  electronic  exchange 
symmetry.  Consequently*  the  total  degeneracy  is  given  by  g*-3(6+(-iyX2/H)*  The  symmetries  are  reversed  for  the  B 
state  as  the  electronic  exchange  symmetry  is  anti -symmetric.  In  this  case  the  total  degeneracy  is  given  by  g*“3(6-(- 
lj'Ml'+i). 

Rotational  energy  transfer  rate  constants  were  derived  from  the  results  of  Lawrence  el  al.t9  Within  the  ground  state* 
the  room  temperature  rate  constants  for  loss  of  X‘  in  the  presence  of  He  and  t)2  were  kmt.x- m*  =  4.7- 10'10  cm3  molecule1 
s'1  and  koulJ(.Qi-  6.3- 1  O' 10  cm3  molecule1  s*\  respectively.  Similarly  for  the  B'  state,  g  Nr  ~  5*5- KT10  cm3  molecule’ * 

s'1  and  3.8  10'10  cm3  molecule'1  s’**  Each  of  these  rate  constants  was  scaled  by  (T/295)172  as  discussed  above. 

The  rate  constants  for  the  transfer  of  population  into  the  levels  X '  and  B  '  from  the  surrounding  rotational  manifolds  were 
then  determined  using  the  principle  of  detailed  balance.  For  the  X  state  the  ratio  was  defined  by  the  expression. 
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In  these  equations  qnit  is  the  rotational  partition  function*  J  is  the  rotational  quantum  number  for  the  X'  or  B'  level,  Bmt 
is  the  rotational  energy  constant  and  is  0*0373679  cm’1  (7.42294-  10'2S  J)  for  the  X  state  and  0.0289937  cm'1 
(5*75944- HP5  J)  for  the  B  state.  The  rotational  partition  functions  for  the  X  and  B  stales,  including  nuclear  spin 
weighting*  are  given  by: 
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Truncation  ofthe  summations  at  7=200  was  found  to  be  a  conservative  approximation  for  temperatures  below  300  K, 

The  absorption  cross  section  <rab![  was  calculated  here  as  opposed  to  the  use  of  measured  values  as  the  experiments 
determining  those  values  were  performed  at  atmospheric  conditions  whereas  the  pressures  used  here  were  1  to  2  orders 
of  magnitude  smaller.  The  absorption  cross  section  is  related  to  the  rate  of  excitation  from  the  X '  stale*  as  given  by  the 
expression: 
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where  g(v)  is  the  Voigt  lineshape  function  given  by; 


Avd  and  Avt  are  the  Doppler  and  Lorentzian  widths;  vand  it  are  the  frequency  and  line  center  frequency  respectively;  vt 
is  the  molecular  velocity  of  the  molecule  in  the  (h  Cartesian  coordinate  direction;  c  is  the  speed  of  light;  Vi  is  the  bulk 
velocity  component  in  the  fh  Cartesian  coordinate  direction;  A  is  the  wavelength  of  the  given  transition;  T  is  the 
temperature;  R  is  the  gas  constant;  m!2  is  the  molecular  mass  of  I2;  T„f  is  the  reference  temperature  for  the  pressure 
broadening  coefficients;  P  is  the  pressure;  and  a,  and  x*  are  the  pressure  broadening  coefficient  and  mole  fraction  of  the 
ih  species  upon  the  given  transition.  Following  the  work  of  McDaniel,  the  pressure  broadening  is  taken  to  be  12 
MHzTTorr  for  the  bath  at  room  temperature,  Tnf  -  298.15  K  and  this  value  is  scaled  by  Tl/2.  The  laser  pulse  properties 
that  determine  p(vj)  are  known  quantities  measured  during  calibration  of  the  laser  source.  The  shape  of  the  pulse  is 
assumed  to  be  Gaussian  in  both  temporal  and  frequency  space.  The  bandwidth  of  the  pulse  was  measured  to  be  22.53 
GHz  and  the  energy  of  the  pulse  was  25- 10‘ 3  J  spread  over  30  ns.  The  Einstein  B  coefficient  for  the  transition  in 
question,  Bv  v%  is  calculated  using: 
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Where  the  squares  of  the  vibrational  overlap  integrals,  |(^(v,)j^(v,,))J  are  the  Franck-Condon  factors  and  the 


transition  dipole  moment  has  been  factored  out  using  the  /f-centroid  approximation.  The  Franck-Condon  factors  used 
here  came  from  a  database  updated  relative  to  the  frequently  cited  work  of  Tellingbuisen15;  the  values  for  the  particular 
transitions  needed  here  are  listed  in  Table  !.  The  values  for  the  rates  of  radiative  decay  for  the  vibrational  manifolds 
within  the  B->X  transition  needed  here  for  determination  the  Einstein  B  coefficient  were  taken  from  the  measurements 
of  Broyer  et  al16  and  Vigue  et  a!17  and  are  listed  in  Table  2. 

The  ordinary  differential  equations  [1J  for  the  loss  and  production  of  the  Xf  X\  B ;  and  B’  states  in  the  presence  of 
the  laser  pulse  are  solved  numerically  using  the  DIVPAG  solution  routine  from  the  IMSL  library  for  each  of  the 
transitions  listed  in  Table  1  within  each  cell  of  the  CFD  dataset.  The  resulting  profile  for  the  temporal  development  of 
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the  B  and  B  '  states  for  each  of  the  pumped  transitions  can  be  used  to  determine  the  number  of  photons  produced  within 
each  by  integrating  the  product  of  the  rate  of  radiative  relaxation,  and  the  total  B  state  number  density,  over  ihe 
period  of  the  camera  gate: 


[18) 

The  total  production  of  photons  traveling  within  the  solid  angle  defined  by  the  light  collection  angle  between  the 
camera  and  the  laser-illuminated  sheet  within  the  gas  defines  the  predicted  collection  of  photons  by  the  CCD  during  the 
camera  gate. 


Results. 

The  3-D  Navier-Stokes  simulation  for  the  PLIF  experiment  nozzle  hardware  shown  in  Figs.  1  and  2  using  the 
computational  grid  shown  in  Fig.  3  was  executed  to  develop  a  prediction  of  the  spatial,  and  if  necessary,  temporal 
variation  of  the  fiowfield.  The  model  was  executed  with  a  timestep  of  l-KT*  sec  for  100,000  time  steps  on  the  high 
resolution  grid  of  51.2  million  grid  cells  after  executing  on  lower  resolution  grids  using  l/8*h  and  1/64*  the  number  of 
grid  cells.  During  the  course  of  the  execution  of  the  simulation  on  each  of  these  varying  resolution  grids,  no  evidence  of 
flow  fluctuation  was  encountered,  and  steady-state  convergence  of  at  least  2  orders  of  magnitude  was  obtained  in  each 
case.  Fig.  5  illustrates  the  predicted  structure  of  the  jet  as  shown  in  an  isosurface  of  constant  12  mole  fraction.  Visible 
within  this  plot  is  the  penetration  and  turn-over  of  the  jet,  as  well  as  the  outline  of  the  counter-rotating  vortex  pair 
embedded  within  the  jet  structure.  Fig.  6  shows  the  temperature  distribution  within  a  2-D  planar  streamwise  cut  through 
the  center  of  the  jet.  Evident  within  this  perspective  is  the  higher  stagnation  temperature  of  the  secondary  fluid  with 
respect  to  the  primary,  400  K  versus  300  K,  shown  by  the  higher  static  temperature  of  the  jet  fluid  for  most  of  the 
duration  of  its’  transit  throughout  the  nozzle.  These  temperature  variations  are  expected  to  impact  the  predicted 
variations  in  photon  yield  from  the  model  for  the  fluorescence  process. 

The  model  for  the  l2  fluorescence  processes,  as  mentioned  above,  can  be  executed  in  a  de-couplcd  manner  from  the 
CFD  model  due  to  the  disparity  in  timescales  between  the  two.  Thus,  the  l2  fluorescence  model  can  be  executed  as  a 
post-processor  for  the  CFD  dataset,  Since  the  planar  L1F  sheets  are  much  smaller  than  the  total  physical  domain  of  the 
nozzle  fiowfield,  a  much  smaller  computational  domain  residing  only  within  the  PLIF  sheet  can  be  created  as  a  subset  of 
the  larger  domain  of  the  CFD  simulation.  To  accomplish  this,  grids  I  mm  thick  matching  the  overall  size  of  the  laser 
sheets  were  created  and  imported  into  the  GASP  software  to  interpolate  from  the  CFD  solution  on  the  original  grid  to 
these  ‘  1  mm  sheet’  grids.  After  interpolation,  the  sheet  grids  with  the  quantities  required  for  input  to  the  l2  fluorescence 
model  including  species  densities  and  mole  fractions,  temperature,  pressure,  and  cell  volumes  were  exported  and  saved 
as  separate  datasets.  The  l2  fluorescence  model  was  then  executed  on  these  sheet  datasets.  The  600  /m  thickness  of  the 
laser  light  is  modeled  within  the  equations  for  the  photon  flux  by  treating  the  profile  as  a  spatial  Gaussian  profile  with  a 
full-width  half-maximum  of  600  pm.  Fig.  7  illustrates  the  location  of  one  of  these  grids  used  for  illuminating  a 
streamwise  cut  norma!  to  the  nozzle  top  wall  through  the  center  of  the  I2/He  injector  orifice. 

The  basic  set  of  principals  underlying  the  model  for  the  12  fluorescence  process  are  that  a  local  ground  slate 
population  of  J2  is  pumped  to  an  excited  state  within  the  B  electronic  state  manifold  and  through  rotational  relaxalion  the 
ground  and  excited  state  populations  are  enlarged  and  diminished  during  the  course  of  the  laser  pumping  process.  Figs. 
8  and  9  demonstrate  the  model  prediction  for  the  temporal  response  of  the  ground  and  excited  states,  X,  X'.  B.  and  B'  to 
the  30  ns  laser  pulse  for  a  location  within  the  He/I2  jet.  Fig.  8  shows  the  temporal  variation  of  the  B'  variable 
representing  ihe  v'.  J  excited  state,  for  each  of  the  lines  pumped.  Transitions  7,  2, 4,  and  6  are  the  strongest  in  terms  of 
population,  respectively,  with  the  remaining  lines  being  significantly  weaker.  The  strength  of  these  lines  over  the 
remaining  lines  is  due  to  a  combination  of  higher  ground  state  populations  and  larger  Einstein  B  coefficients.  Fig.  9 
shows  the  temporal  variation  of  the  X,  X',  B,  and  5'  variables  for  Line  7  over  the  course  of  the  laser  pulse.  The  effect  of 
the  filling  of  the  ground  rotational  state  by  the  adjacent  rotational  states  is  apparent  here  as  the  X'  population  decreases 
by  only  10%  during  the  course  of  the  pulse.  However,  the  total  population  that  is  generated  in  the  B  and  £'  variables  for 
the  excited  state  is  20  times  larger  at  the  peak  than  that  lost  from  the  X\  thus  clearly  illustrating  the  importance  of  the 
rotational  energy  transfer  from  the  adjacent  ground  states  in  the  generation  of  the  excited  states  and  in  turn  photons  from 
fluorescence. 

The  radiative  decay  rate  determines  the  photon  flux  from  the  excited  states  in  the  model.  The  total  fluorescence 
decay  lifetime  (reciprocal  of  the  decay  rate  defined  by  Eq.  [3]),  is  measured  in  the  PLIF  experiment,  providing  an  avenue 
for  validation.  Fig.  10  shows  a  comparison  between  a  plot  of  the  lifetime  as  measured  in  the  PLIF  experiment  with  that 
predicted  by  the  fluorescence  model  within  the  PLIF  sheet  illustrated  in  Fig.  7.  The  lifetime  values  within  the  He/1 2  jet 
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were  on  the  order  of  50  ns  with  lengthening  of  the  lifetime  to  order  IO0  ns  further  downstream  as  the  collisional 
influence  decreased  due  to  decreasing  temperature  and  pressure.  Overall,  the  comparison  between  the  experiment  and 
model  lifetime  values  is  good,  validating  the  quantitative  and  quantitative  aspects  of  the  lifetime  calculation  portion  of 
the  model. 

Direct  comparisons  between  PLIF  images  from  the  experiments  and  images  generated  using  the  I2  fluorescence 
model  are  shown  in  Figs.  11 -15.  The  model  images  were  generated  by  summing  the  total  photon  yield  in  the  direction 
orthogonal  to  the  sheet,  and  dividing  by  the  cross-sectional  area  of  the  cell  faces  parallel  to  the  sheet  to  convert  absolute 
photon  counts  into  photon  fluxes  or  effectively  intensity.  Fig.  11  shows  the  comparison  within  the  sheet  through  the 
injector  orifice  parallel  to  the  flow  direction.  The  primary  jet  is  clearly  illuminated  with  the  highest  intensities  in  the 
image  with  a  secondary  structure  associated  with  the  counter-rotating  vortex  pair  (CRVP)  seen  between  the  primary  path 
of  the  jet  fluid  and  the  nozzle  wall.  In  the  model  generated  image,  the  region  between  the  primary  jet  path  and  the 
structure  associated  with  the  CRVP  is  not  illuminated  as  brightly  as  in  the  experiment  image.  This  difference  may  be 
due  to  an  inaccuracy  of  the  model,  or  may  be  due  to  saturation  effects  within  the  experiment  image  that  are  not  captured 
here.  In  all  images  shown  here,  the  lightest  or  whitest  grayscale  shades  are  associated  with  the  brightest  fluorescence, 
and  he  maximum  in  the  contouring  scale  is  set  to  this  value.  However,  in  the  experiment  images,  some  pixels  may  have 
been  saturated  during  the  course  of  the  light  collection,  and  would  not  accurately  register  the  highest  fluorescence 
intensities,  thus  shifting  the  ratio  of  highest  to  lowest  intensity  recorded  downward  and  shifting  the  grayscale  contouring 
in  the  image.  The  end  effect  would  be  a  brightening  of  lower  intensity  regions  relative  an  image  without  saturation.  As 
the  saturation  characteristics  of  the  CCD  camera  used  to  collect  the  experiment  images  are  not  built  into  the  I2 
fluorescence  model,  differences  between  the  model  and  experiment  images  with  some  saturation  will  be  found. 
Nonetheless,  the  comparison  between  the  jet  penetration  and  trajectories  is  very  good  between  both  images. 

Figs.  12-15  show  images  from  light  collection  in  which  the  planar  laser  sheet  interrogates  the  fiowfield  at  a  45° 
angle  with  respect  to  the  primary  flow  direction  as  illustrated  in  Fig.  L  Figs.  12  and  13  show  comparisons  between  the 
experiment  images  at  the  trailing  edge  of  the  injector  orifice  with  two  separate  images  from  the  model.  The  first 
comparison  in  Fig.  12  uses  a  fluorescence  model  image  in  which  the  predicted  number  of  photons  generated  during  the 
laser  pulse  is  integrally  summed  in  the  direction  orthogonal  to  the  sheet  to  mimic  the  effect  of  the  finite  sheet  thickness 
of  600  /jiti.  While  both  images  show  clear  illumination  of  the  edges  of  the  fluid  within  CRVP,  the  model  shows  more 
illumination  within  the  center  region  of  the  jet  fluid.  This  difference  could  be  due  to  a  thinner  sheet  of  fluid  being 
illuminated  within  the  experiment  than  the  indicated  sheet  thickness  of  600  /mt,  an  important  consideration  when 
imaging  structures  having  the  same  scale  as  the  sheet  thickness.  Fig.  13  compares  an  image  from  the  model  using  only  a 
single  plane  at  the  center  of  the  sheet  instead  of  the  integrally  summed  technique.  The  model  prediction  shows  sharp 
resolution  of  the  CRVP  structure  as  do  the  original  and  the  experiment  images,  hut  more  fluid  remains  resolved  within 
the  CRVP  than  in  the  experiment.  Thus,  it  is  possible  that  slight  differences  in  the  location  of  the  sheet  between  the 
experiment  and  the  model  could  be  contributing  to  the  illumination  of  slightly  different  segments  of  the  same  fluid 
structure.  Figs.  14  and  15  show  comparisons  further  downstream,  at  locations  3  and  7  cm  downstream  of  the  injector 
orifice.  The  model  images  used  the  integral  summing  of  the  photon  contributions  through  the  thickness  of  the  sheet,  and 
very  good  agreement  is  demonstrated  between  the  experiment  and  model  images,  dearly  showing  nearly  identical 
structure  of  the  CRVP  and  associated  fluid  within  the  jet.  The  model  images  are  slightly  darker  than  the  experiment 
images  as  the  grayscale  contouring  is  fixed  to  that  used  in  the  injector  imaging  in  Fig.  12,  and  tends  to  confirm  the 
notion  regarding  the  shift  in  the  grayscale  contouring  in  the  experiment  images  due  to  saturation. 

Summary  and  Conclusions. 

In  conclusion,  the  results  show  here  demonstrated  the  development  and  application  of  a  model  for  the  fluorescence 
of  1 2  in  the  presence  of  a  pulsed  laser  source.  Very  good  agreement  between  the  model  and  the  experiment  PLIF  images 
was  found,  but  at  the  same  time  illustrated  the  challenges  for  comparing  high  resolution  imagery  and  high  resolution 
CFD  simulations. 

Future  work  will  focus  on  detailed  examinations  of  the  quantification  of  CCD  pixel  saturation,  precise  quantification 
of  the  laser  sheet  positioning  and  width,  and  parameterizations  of  mixing  studies. 
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Tabic  1 .  Data  for  the  transitions  used  in  the  i2  fluorescence  model 
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Table  2.  Lifetime  data  used  in  the  12  fluorescence  model. 
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Figure  1,  Representation  of  the  PLIF  experiment  configuration  simulated  here.  The  laser  sheet  is  passing 
through  the  experiment  at  a  45®  angle  and  the  camera  is  oriented  orthogonally  to  the  illuminated  sheet. 


Figure  3.  Computational  grid  used  in  3-D  Navier-Stokes  simulation  of  the  experiment  hardware  fiowfield 
in  Fig.  I,  Grid  resolution  shown  ts  1/5 1 2th  of  that  used  in  the  actual  simulation. 
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Figure  2.  Detailed  drawing  for  the  l2  injector  orifice  used  in  the  PLIF  nozzle  experiment  hardware. 


for  the  nozzle  simulation. 


21 


,tk 

f  '  X 


Figure  5*  Isosurface  plot  for  constant  mole  fraction  of  12t  value  =4*1  Q'3,  from  the  3-D  Navier-Stokes 
simulation  of  the  PLIF  nozzle  experiment  flowfield.  The  outline  of  the  counter-rotating  vortex  pair  is 
visible  in  the  upper  surface  structure  downstream  of  the  orifice. 
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Figure  7.  Illustration  of ‘sheet1  sub-domain  used  for  execution  of  the  12  fluorescence  model,  in  blue,  with 
respect  to  the  larger  computational  domain  used  in  the  3-D  Navier-Stokes  simulation.  Values  are 
interpolated  from  the  Navier-Stokes  simulation  domain  to  the  ‘sheet1  domain. 
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Figure  6*  Temperature  (K)  contours  within  a  2-D  planar  cut  through  the  center  of  the  jet  within  the  3-D 
Navier-Stokes  simulation  of  the  PUF  nozzle  experiment  flowfield.  Note  that  the  jet  fluid  remains  warmer 
than  the  surrounding  fluid. 
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Figure  8.  Predicted  temporal  variation  of  the  B'  states  for  a  data  point  within  the  l2  jet  illustrating  the 
variation  in  excitation  during  the  laser  pump  from  line  to  line. 
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Figure  9,  Predicted  temporal  variation  of  Line  7,  i/-l,  7' -75,  J  =  74  from  the  l2  fluorescence 

model.  Note  that  the  Y  axis  for  the  X  states  variable  is  to  the  right  and  that  for  the  remaining  variables  is  to 
the  left. 


Figure  11.  Comparison  of  experiment  PDF  image  with  fluorescence  model  generated  image  within  the 
streamwise  laser  sheet  shown  in  Fig.  7 
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Figure  10.  Comparison  of  the  fluorescence  lifetime  as  measured  in  the  PLIF  nozzle  experiment,  top  image, 
with  (hat  from  the  fluorescence  model  Note  that  the  experiment  includes  reflection  of  the  fluorescence 
light  off  of  the  nozzle  wall. 


Figure  12+  Comparison  of  experiment  PLIF  image  with  fluorescence  model  generated  image  within  a  sheet 
interrogating  the  flowfldd  at  a  45*  angle  with  respect  to  the  flow  axis  at  the  trailing  edge  of  the  injector 
orifice.  Note  that  the  model  generated  image  is  integrally  averaged  across  the  width  of  the  sheet. 
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Figure  13,  Comparison  of  experiment  PL1F  image  with  fluorescence  model  generated  image  within  a  sheet 
interrogating  the  flowfield  at  a  45°  angle  with  respect  to  the  flow  axis  at  the  trailing  edge  of  the  injector 
orifice.  The  model  image  is  extracted  from  the  plane  at  the  center  of  the  sheet,  and  does  not  include  light 
contributions  from  the  full  thickness  of  the  sheet. 


Figure  15.  Comparison  of  experiment  PLIF  image  with  fluorescence  model  generated  image  within  a  sheet 
interrogating  the  flowfield  at  a  45°  angle  with  respect  to  the  flow  axis  7  cm  downstream  of  the  injector 
orifice. 
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Figure  14.  Comparison  of  experiment  PLIF  image  with  fluorescence  model  generated  image  within  a  sheet 
interrogating  the  flowfield  at  a  45°  angle  with  respect  to  the  flow  axis  3  cm  downstream  of  the  injector 
orifice. 
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